function [ equal ] = isEqualRel( x, y, error )
% Compara dos matrices, usando el error
    relTol = error;   % This one too!
    absError = x(:)-y(:)
    relError = absError./x(:)
    relError(~isfinite(relError)) = 0  % Sets Inf and NaN to 0
    equal = all( abs(relError) < relTol )

end

